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Abstract 

The triple-differential dijet cross-section, d 3 a ep /dQ 2 dE t dxif s , is measured with the HI 
detector at HERA as a function of the photon virtuality Q 2 , the fraction of the photon's 
momentum carried by the parton entering the hard scattering, xiy ets , and the square of the 

2 

mean transverse energy, E t , of the two highest E t jets. Jets are found using a longitudinal 
boost-invariant kx clustering algorithm in the 7*p center of mass frame. The measurements 
cover the ranges 1.6 < Q 2 < 80 GeV 2 in virtuality and 0.1 < y < 0.7 in inelasticity y. 
The results are well described by leading order QCD models which include the effects of 
a resolved component to the virtual photon. Models which treat the photon as point-like 
fail to describe the data. An effective leading order parton density for the virtual photon is 
extracted as a function of the photon virtuality, the probing scale and the parton momentum 
fraction. The x 7 and probing scale dependences of the parton density show characteristic 
features of photon structure, and a suppression of this structure with increasing Q 2 is seen. 
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1 Introduction and Motivation 



It is well-established that the real photon has a partonic structure both through measurements 
in two-photon collisions at electron-positron colliders [[IJ] and through measurements of the 
photoproduction of jets at HERA |3|, |J]. These data have been used to determine universal 
parton densities for real photons. The dynamical evolution of the partonic structure of the 
photon as it becomes virtual is described in a QCD framework both for deep-inelastic positron- 
proton (ep) collisions and for high transverse momentum (P t ) processes in two-photon and 
photoproduction reactions [|5|] . Experimental data with target photons of sizeable virtuality in 
two-photon collisions are, however, sparse [§]. The ep collision data at HERA on the other 
hand are available over a wide range of photon virtualities, Q 2 , from photoproduction to high 
Q 2 deep-inelastic scattering (DIS) and are sensitive to any photon structure [|7|, || |9|]. 

The production of high transverse energy (E t ) jets in ep collisions is dominated by processes 
in which a single space-like photon carries the momentum transfer from the incident positron. 
When this photon is quasi-real, i.e. in photoproduction processes, two types of interaction can 
be distinguished in leading order: direct processes in which the photon couples as a point-like 
object to a parton out of the proton and resolved processes in which it develops a partonic struc- 
ture prior to the collision. In the latter case, a parton out of the photon, carrying only a fraction 
x 7 of the photon momentum, enters the hard scattering process leading to the production of 
jets. Examples of these two types of process are shown in figure 1. The photoproduction jet 
cross-section is therefore sensitive to the density of partons in the photon and the latter can be 
measured. A natural choice of the scale, P t , at which this photon structure is probed is given by 
the E t of the jets with respect to the axis in the 7^ centre of mass system (cms). 

The diagrams in figure 1 are equally applicable to processes where the exchanged photon 
is highly virtual. In such deep-inelastic scattering processes, the production of high E t jets is 
usually dominated by direct processes. However, as long as the photon is probed with suffi- 
ciently high resolution, i.e. if P 2 ^> Q 2 (with P t again defined by the E t of the jets in the 
7*p frame), it may still be possible to have interactions in which the cross-section factorises 
and the photon structure is resolved [|7], 0, [ITJ [TTJ, [T^p. It is then a prediction of perturbative 
QCD that the parton densities of virtual photons become suppressed as Q 2 increases at fixed 
P 2 . The partonic structure in the photon becomes simpler, until only the direct coupling to a qq 
pair remains [0, |9|, |10|, [n], [T^]. The concept of virtual photon structure thus provides a unified 



description of high E t jet production over the whole range of Q 2 corresponding to a smooth 
transition between photoproduction and DIS processes. It is implicit in this concept that the 
structure is universal. Future comparisons with, for example, 7*7 data may establish whether 
this is the case. 

In a previous publication JT3|], it was shown that the single inclusive jet cross-section in low 
Q 2 deep-inelastic scattering can indeed be described by models which include a contribution 
from resolved virtual photons which is suppressed with increasing Q 2 . 

The dijet cross-section has also been measured in low Q 2 deep-inelastic scattering processes 
and has been found to be best described with predictions in which the effects of virtual photon 
structure are included JT3|] . 
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With dijet events, it becomes possible to reconstruct the variable x 7 and hence extract an 
effective parton density for the photons. Such a measurement has recently been made using 
photoproduction events [t3]] . 

In this paper we extend the latter studies and investigate the evolution of the effective parton 
density with Q 2 as well as the P 2 of the partons. We begin by measuring the triple-differential 
jet cross-section, d 3 a ep /dQ 2 dE t dx^, where E t is the mean of the transverse energies of the 
two highest E t jets measured in the YP centre of mass frame, and x^ ts is the value of x 7 as 
estimated from the jets. Jets were found using the inclusive algorithm [|T6|]. The measured 
cross-section is compared to simulations with LO matrix elements which include models for 
the virtual photon structure. 

The cross-section is then used to extract a leading order effective parton density as a function 
of x 7 , Q 2 and P 2 . The observed shape, scaling behaviour and virtuality dependences of this 
effective parton density are compared with various parameterisations based on predictions from 
perturbative QCD. 



2 The HI Detector 

The HI detector is described in detail in []17[]. In this analysis, we make particular use of 



the SPACAL Q180 calorimeter for detection and identification of the scattered positron. The 



hadronic energy flow is measured with the Liquid Argon (LAr) [TlPp and SPACAL calorime 



ters. The central and forward tracking detectors are used to reconstruct the event vertex and to 
supplement the measurement of hadronic energy flow made with the LAr and SPACAL. 

We use a coordinate system in which the nominal interaction point is at the origin and the 
incident proton beam defines the +z direction. The polar angle 9 is defined with respect to the 
proton direction. 

The central tracking system consists of two concentric cylindrical drift chambers, coaxial 
with the beam-line and centered about the nominal interaction vertex. Its polar angle cover- 
age, 15° < 9 < 165°, is complemented by that of the forward tracker, 7° < 9 < 25°. The 
central tracker is interleaved with drift chambers providing measurement of the z coordinates 
of the tracks. The tracking detectors are immersed in a 1.15 T magnetic field generated by a 
superconducting solenoid which surrounds the LAr. The LAr is a finely grained calorimeter 
covering the range in polar angle 4° < 9 < 154° with full azimuthal acceptance. It con- 
sists of an electromagnetic section with lead absorbers, 20-30 radiation lengths in depth, and a 
hadronic section with steel absorbers. The total depth of the calorimeter varies between 4.5 and 
8 hadronic interaction lengths. The energy resolution is a(E)/E m 0.12/V^©1% for positrons 
and a(E)/E « 0.5/ \f~E~ © 2% for hadrons (E in GeV), as measured in test beams tf2X$\. The 
absolute energy scale is known to a precision of 3% for positrons and 4% for hadrons. 

The SPACAL is a lead/scintillating-fibre calorimeter which covers the angular region 
153° < 9 < 177.8°. It contains an electromagnetic and a hadronic section. The former has 
an energy resolution of 7.5%/ \fE © 2.5%. The energy resolution for hadrons is ~ 30%/ \/E. 
The SPACAL provides the main trigger for the events (see section ^) in this analysis. The tim- 
ing resolution of better than 1 ns in both sections of the SPACAL is exploited in the trigger to 
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reduce proton beam induced background. The energy scale uncertainty of the electromagnetic 
section is 2% and that of the hadronic section is 7%. The backward drift chamber (BDC) system 
in front of the SPACAL spanning the region 151° < 9 < 177.5° provides track segment infor- 
mation to improve positron identification in the SPACAL. In conjunction with the event vertex 
determination from the central and forward track detectors, it gives a precision measurement of 
the positron scattering angle of 1 mrad. 

The luminosity determination is based on measurement of the ep — > ep7 Bethe-Heitler 
process. The positron and photon are detected in the electron tagger located at z = —33.4 m 
and photon tagger at z = — 103 m, respectively. Both consist of crystal Cherenkov calorimeters 
with a resolution of cr(E) / E rs 0.1/ \[E. The integrated luminosity was measured to a precision 
of better than 2%. 



3 Theoretical Models and Simulations 



3.1 Monte Carlo Models 

The analysis uses simulated events both to correct the measured cross-sections for detector 
effects and in order to compare the data with the predictions of various theoretical models. The 
various combinations of Monte Carlo simulations and parton densities used in this analysis are 
summarised in table [I]. 

The HERWIG flnfl and RAPGAP [[22]] Monte Carlo models are both able to simulate the 
direct and resolved production of dijets by virtual photons. In both models the hard scattering 
process is simulated in leading order (LO), regulated with a minimum P t cut-off, P™ n , and 
supplemented by initial and final- state parton showers. In the simulation of resolved processes, 
the equivalent photon approximation is used for the flux of transversely polarised photons and 
on-shell 2^2 matrix elements are taken. The longitudinal flux is not included. Exact eq — > 
eqg and eg — ► eqq matrix elements are used for the direct photon component. 

HERWIG has parton showering based on colour coherence and uses the cluster model for 
hadronisation [|23|]. Some tuning of the parton showering is possible and we use two different 
settings for the scaleQ at which the parton showering is terminated. In addition to the hard 
scattering process, HERWIG can also model the additional soft underlying activity in the event 
(SUE) which is necessary to describe the observed energy flow in and around the jets [|]]. Such 
soft particles (uncorrected with the hard scattering process) may be produced via soft remnant- 
remnant interactions. The probability that a resolved event contains soft underlying activity 
was adjusted in the simulation. Event samples with 0% and 100% SUE were mixed to simulate 
different probabilities. No SUE was introduced into the direct sample. 



RAPGAP uses a leading-log parton shower approach and the LUND string model Q24J ] for 
hadronisation. It contains no mechanism for simulating additional soft underlying activity in 
the events. 



'Specified by the parameter QSPAC. 
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The simulations have interfaces to a variety of parameterisations of photon and proton parton 
density functions (PDF). The factorisation scales of the proton and photon parton densities were 
both set equal to the P t of the scattered partons with respect to the 7^ axis in the cms. GRV 



parton densities [g5j] were used for the proton. When correcting the data for detector effects, 



higher order (HO) versions of the parton densities and the 2-loop expression for as were used. 



As has also been noted by the ZEUS collaboration Q260, with this configuration it is necessary 



to re-scale the HERWIG predictions by a factor of 1 .7 in order to describe the data. 

Several models for the virtual photon parton densities were considered. The Drees-Godbole 
model (DG) [ fTT| , [12[ ] starts with real photon parton densities [[27]] and suppresses them by a 



factor £ which depends on Q 2 , P 2 and a free parameter, u, which controls the onset of the 
suppression: 



In 



C(Q 2 ,P?^) = ^±^ (1) 

r t -+- U) 



U 2 



Quark densities in the real photon are suppressed by £ and the gluon densities by C 2 . This 
ansatz, based on the analysis in [jT2p, is designed to interpolate smoothly between the leading- 
logarithmic part of the real photon parton densities, ~ \n(P 2 / Aq CD ), and the asymptotic do- 
main, P 2 ^> Q 2 ^> A.q CD , where the photon density functions are predicted by perturbative 
QCD to behave as ~ ln(P 2 /Q 2 ). In this model, the shape of the x 7 distribution evolves with 
Q 2 because of the different suppression of the quark and gluon densities. 



In the models of Schuler and Sjostrand (SAS) []28[], the virtual photon parton densities are 
decomposed into a non-perturbative component modelled by vector meson dominance (VMD) 
and a perturbative anomalous component. As Q 2 increases, the VMD component is rapidly sup- 
pressed, ~ [m v / {rn v + Q 2 )] 2 , whereas the anomalous part has a slower logarithmic suppression 
and is again designed to approach the exact QCD predictions in the P 2 Q 2 Aq CD region. 
There are four models, SAS-1M, SAS-2M, SAS- ID and SAS-2D which differ in their choice 
of factorisation scheme (DIS (D) or MS (M)) and the scale at which the evolution is started (0.6 
or 2.0 GeV indicated by the 1 or 2 in the name, respectively). 

The events used for the correction of detector effects were processed through a full simu- 
lation of the HI detector. The HERWIG event sample which we use as our main model for 
the corrections, contains approximately 3 times as many events as the selected data sample. 
The statistics of the RAPGAP sample, which we use in the estimation of systematic errors, are 
comparable to that of the data. 



4 Event Selection 



The analysis is based on positron-proton collision data collected by the HI detector at HERA 
in 1996 and corresponds to an integrated luminosity of ~ 6 pb~ . During this period, 820 GeV 
protons collided with 27.5 GeV positrons. The events were triggered by an energy deposition 
exceeding a threshold energy in the electromagnetic section of the SPACAL of ^2.5 GeV at 



7 



Model name 


a s 


Proton PDF 


7* PDF 


TDmin 

r t 

(GeV) 


SUE 


QSPAC 
(GeV) 


HERWIG(HO)/DG 


2-loop 


GRV-HO 


GRV-HO*DG 


3 


10% 


1.0 


(xl.7) 






(uj = 0.2 GeV) 








RAPGAP(HO)/SAS-2D 


1-loop 


GRV-HO 


SAS-2D 


3 


- 


- 


HERWIG(LO)/DG 


1-loop 


GRV-LO 


GRV-LO*DG 

(u = 0.2 GeV) 


2 


5% 


2.0 


RAPGAP(LO)/DG 


1-loop 


GRV-LO 


GRV-LO*DG 

(u = 0.2 GeV) 


2 






RAPGAP(LO)/SAS-2D 


1-loop 


GRV-LO 


SAS-2D 


2 






RAPGAP(L0)/SAS-1D 


1-loop 


GRV-LO 


SAS-1D 


2 







Table 1: Description of models. The first two models are used to correct for detector effects. 
Note that the predictions from the first model were scaled by a factor of 1.7. The remainder are 
used in the unfolding of the diparton cross-section from the corrected jet cross-section and for 
comparison with the corrected cross-section and effective parton density. 



large radii or 5.7 GeV at large radii, provided this was accompanied by at least one track in the 
central tracking device with transverse momentum P t ^0.8 GeV and a well-defined interaction 
vertex. The efficiency of this trigger is typically ~ 90% and has been measured over the relevant 
kinematic range to a precision of 5%. 

Additional criteria were then applied to reduce background events and to ensure that the 
events were well-measured. The scattered positron was identified as follows. An energy cluster 
in the SPACAL was required to have a radius of less than 3.5 cm, consistent with being produced 
by a positron. In order that it be well-contained, the cluster was required to be at a radius greater 
than 8 cm, to have less than 10% of its energy in the innermost cells of the calorimeter and to 
have < 0.5 GeV of energy deposited in the cells of the hadronic section immediately behind 
the cluster. Finally, a track in the BDC was required such that the radii of the SPACAL cluster 
and BDC track differed by no more than 3.0 cm. 

The reconstructed event vertex was required to lie within 30 cm of its nominal location 
in z. In order to further reduce photoproduction background, we required that ^^E 1 — P%)> 
which is expected to approximately equal twice the electron beam energy, E e , for DIS events, 
lies in the range 40-65 GeV. The sum is taken over all calorimeter clusters supplemented by 
tracking information. This procedure corrects for energy loss in the passive material in front of 
the calorimeter. Each track was allowed to contribute a maximum of 0.35 GeV to avoid double 
counting with the calorimetric energy. 

The inelasticity, calculated from the energy, E' e , and polar angle, 9', of the scattered positron: 

y e = l-E'jE e sm 2 {e'/2), (2) 

was restricted to the range 0.1 < y e < 0.7. The upper limit corresponds to a minimum scattered 
positron energy of ~ 8 GeV. Small y e values, where the measurement is limited by the energy 
resolution of SPACAL, were excluded. The positron energy and scattering angle were also used 
to calculate the virtuality, Q 2 , of the photon: 

Q 2 = AE' e E e cos 2 (9'/ 2). (3) 
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The inclusive k t clustering algorithm [[16]] was applied to find jets in the j*p cms frame with 



the boost defined by the scattered positron's energy and scattering angle. Tracking information 
was used to improve the reconstruction of the E t of the jets as described above. The calorimeter 
energy clusters were treated as massless objects and the tracks were assigned the pion mass. 
The clustering of these final state objects into jets uses as the distance measures: 

d tj = min (PI 2 , Pi 2 )AR (4) 
where AR = ^/{rf - r]i) 2 + (ft - (pi) 2 (5) 
and di = P\ 2 (6) 

between objects % and j and between the i'th object and the beam respectively. The P t \ rf and 
(f) 1 are the transverse momenta, pseudorapidities given by rf = — lntan(# l /2) and azimuthal 
angles of the objects, respectively. At each iteration of the algorithm, the smallest distance 
measure is determined. If this is one of the d^, then the objects i and j are combined into a 
massless object using the P t weighting scheme: 

p t = p: + pi (7) 

i P i + j P j 

= (8) 



Pt 

Whenever a di is the smallest distance measure, the z'th object defines a completed jet and is 
excluded from further iterations. The iterations terminate when all objects have been assigned 
to jets. The reconstructed jets are massless. 

Events were required to contain at least two jets found by the above algorithm. Only the 
two highest E t jets are used in the analysis. Events were accepted if the two highest E t jets in 
the event (called in the following jet 1 and jet 2) satisfied the following criteria: 

\rf etl -rj jet2 \ < 1.0 (10) 
-3.0<r/<-0.5 (11) 
E\ > 30 GeV 2 (12) 

■Mrs ^~pt~ < 0.25 (13) 

{Ef l + Ef 2 ) 

where fj and E t are the mean pseudorapidity and mean transverse energy of the two highest E t 
jets. Et and r] are always given with reference to the 7*p cms with the proton direction defining 
the positive ^-axis. The same cuts, applied to the jets after correction for detector effects, serve 
to define the cross-section. The cuts are the same as those used in reference ||3j]. The first two 
cuts help to ensure that the jets are confined to a region of the detector where they are well- 
measured and therefore that x 7 can be well-determined. The restrictions on the difference in t] 
and Et of the jets reduce the probability of misidentifying a part of the photon or proton remnant 
as one of the high E t jets. The constraints on E t are such that neither jet has E t < 4 GeV and 
the sum of the jet E^s is always >11 GeV. In 10% of the selected events there is a third jet with 
E t > 4 GeV. Although we do not do so here, the asymmetry in the jet selection makes possible 
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comparison with NLO QCD calculations which become unstable for symmetric jet cuts if 
both jets are near their common lower limit. 

After application of these selection criteria we obtained a sample of approximately 12,000 
dijet events with ~Et > 30 GeV 2 spanning the Q 2 range 1.6 < Q 2 < 80 GeV 2 . Diffractive 
events were not explicitly excluded. The residual background from photoproduction is negligi- 
ble. 



5 Measurement of the Triple Differential Dijet Cross-section 

We first study the dependence of the dijet cross-section on the variables xi^ ts , E T and Q 2 . The 
jet-based variable, xi^ ts , is related to the true a; 7 of the events: 

E ( Ejet - ^ ei ) 

= ^ : „ • (14) 

h 

Energies and momenta are measured in the YP cms an d with respect to the j*P axis. The sum 
in the denominator is over all final state particles (except for the scattered positron). As follows 
from the conservation of energy and longitudinal momentum, x^ ts is equal to the true a; 7 for 
leading order dijet production. 

In each event, an estimate of x 3 y ets was made using the reconstructed energies and longi- 
tudinal momenta of the two highest E t jets in the 7*p cms. This quantity is referred to as 
x rec below. sum m the denominator was taken over all reconstructed objects in the event 
(calorimeter clusters supported by tracking information) except for those associated with the 
scattered positron. 



5.1 Correction of the Data for Detector Effects 

The measured dijet cross- sections were corrected for detector acceptance and resolution ef- 
fects in the kinematic domain specified in the previous section. An iterative two-dimensional 
Bayesian unfolding technique [ ^0| ] was applied to distributions of x'^ c and E T rec in separate 
ranges of reconstructed Q 2 . The correlations between the measured variables {x^ ec , E^ rec } 

and the corrected variables {x^ ets , E T } were obtained using events generated by the HER- 
WIG(HO)/DG model (see Table 1). The generated events were subjected to a detailed sim- 
ulation of the HI detector. The x^ ts resolution is approximately 12%, independent of the xif ts , 

2 2 2 

E t and Q 2 values. The resolution in E t is ~ 14% in the highest E t range reducing in pre- 
cision to ^40% in the lowest. The resolution in Q 2 is 8% at low Q 2 and 2% at high Q 2 . A 
bin-by-bin correction was then applied for the Q 2 dependence. The measurement error on Q 2 is 
much smaller than the bin size and migrations are negligible. After unfolding, the correlated er- 
ror between bins in the unfolded distributions was < 40%. The systematic errors are described 
below. Here as elsewhere in this analysis, the various systematic errors are added in quadrature 
and unless otherwise stated, the numbers represent average values for the errors. 
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Model Dependence. The largest sources of systematic error arise from model depen- 
dences, in particular from the choice of parton showering and hadronisation models. 
These were estimated by comparing the results of unfolding using HERWIG with those 
from unfolding with the RAPGAP(HO)/SAS-2D simulation and then by assigning the 
full difference as the error symmetrically to the measurement. The estimate of ~ ±20% 
includes the small contribution from the choice of input parton densities. 

Stability of Unfolding. By varying the number of iterations used, we estimated that un- 
folding instabilities result in a 5% systematic uncertainty in the unfolded cross-sections. 

Absolute Energy Scales. The uncertainties in the LAr and SPACAL calorimeter hadronic 
energy scales lead to 12% and 1% systematic errors in the results, respectively. The 
uncertainty in the SPACAL electromagnetic energy scale yields a 4% systematic effect. 

Trigger Efficiency. The trigger efficiency uncertainty results in a 7% error. 

Radiative Corrections. Radiative corrections have been estimated to result in a change in 
the cross-sections which is typically less than 5%. As no simulation of resolved photon 
processes which includes radiative corrections is available, the estimate is based on direct 
events only. The data are not corrected for this and it is not included in the systematic 
error. 

Soft Underlying Event. The unfolding procedure might also be influenced by the soft 
underlying event. In Figure |2| we show the measured transverse energy flow about the 
two highest E t jets. The flow is calculated from the energy clusters in a strip of ±1 unit 
of i] with respect to the jet axes as a function of the difference between the <p of the clusters 
and the <p of the jet axis. Tracks are also included, modified as described in section 4. The 
distribution is made for each of the two highest E t jets in the events and separated into 
different ranges of jet E t and r]. Because the two jets are constrained to be close together 
in r], the energy flow associated with the other approximately back-to-back jet is clearly 
visible. Note that the flow produced by this other jet is dependent on event topology as 
well as on the jet properties. As the two jets are not precisely back-to-back, by orienting 
the flow in <\> such that the axis of the other jet is always to the left, we expose a wider, 
more clearly visible pedestal region to the right. 

The flow is compared with the HERWIG(HO)/DG simulation with 10% soft-underlying 
event and with predictions from RAPGAP(HO)/SAS-2D which includes no model for the 
soft underlying event. Both simulations give a good description of the energy flow in the 
core of the jets. Neither model is able to describe the energy flow in the pedestal region 
for all ranges of E t and r]. For i] < —2.0, the pedestal is well-described by RAPGAP. 
In the forward region, rj > -1.25, the data lie between the RAPGAP and HERWIG 
predictions. Although HERWIG overestimates the data in this inclusive plot, we find that 
10% soft underlying event is needed to account for the pedestal observed for Q 2 $ 8 GeV 2 
(not shown separately) [[Bp. The unfolding was repeated using HERWIG(HO)/DG with 
0% and 15% soft underlying event. Differences in the resulting distributions were found 
to be ~ 3% and are included in the systematic errors. The jet pedestal has more influence 
on the measurement of the effective parton density and will be discussed again later in 
section 6. 
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5.2 Discussion of the Results 



The corrected triple-differential cross-section for 0.1 < y < 0.7 and for jets satisfying the 
criteria given in equations 10 - 13 is given in Tables [I] and |3| and is shown in Figures ^ |] 
and |5| in various projections. The three Figures show the cross-section as functions of xif ts , 

E t and Q 2 , respectively. In each case, the distributions are shown for ranges of the other two 
variables. The data are depicted as points with error bars. The error bars show the statistical and 
systematic errors added in quadrature. Note that the systematic errors are dominant everywhere. 
The absence of a data point indicates that no measurement was made in that bin because of 
insufficient statistics for the unfolding. 

In Figure [|, the a£ ets distributions can be seen to peak towards a£ ets = 1, where the direct 
photon contribution is expected to be concentrated. There is a strong decrease in the cross- 
section with increasing E t . As Q 2 increases at fixed E t , the cross-section decreases and, for 
E t < 150 GeV 2 , the relative contribution from resolved photons, in the region x^ ts ^0.75, 

can be seen to diminish. For E t 2 > 150 GeV 2 , only the highest x^ ts point has been measured. 
Note that the reduction in the cross-section in the lowest x^ ts bin is a consequence of the cuts 

2 

in 7) and E t . 

The data are compared with predictions of models with LO matrix elements and parton 
densities. The HERWIG(LO)/DG simulation is shown for two choices of the Q 2 suppression 
factor, uj. The predictions for the direct only contribution are also shown. With uj — 0.1 GeV, 
the HERWIG(LO)/DG model gives a reasonable description of the cross-section throughout 
the Q 2 -E t range. Increasing uj to 0.2 GeV leads to an overestimation in the low Q 2 , low E t 
regions. The value of uj which best describes the data, however, depends on the frequency 
of soft underlying events. The RAPGAP(LO)/DG model prediction (not shown in this Figure), 
which contains no soft underlying event, requires uj ~ 0.2 GeV in order to describe the data. As 
Q 2 increases, both HERWIG models tend to underestimate the cross-section for intermediate 
x jets values and overestimate it at high x^ ts The direct only contribution is able to describe the 
data in the highest x^ ts bin but underestimates the data for x^ ts ^0.75 except possibly in the 
highest Q 2 range. 

2 

The E t dependence of the cross-section is shown in Figure |J It is compared with pre- 
dictions from RAPGAP(LO) using three different choices for the photon parton densities. In 
the highest xi, ets range, where direct processes are expected to dominate, all the models give 
a good description of the data. Elsewhere, the models provide a spread of predictions and no 
single one is preferred. The Drees-Godbole model tends to overestimate the cross-section for 
1.6 < Q 2 < 3.5 GeV 2 . For Q 2 > 3.5 GeV 2 and 0.45 < x 7 < 0.75, all three models tend to 
underestimate the data but are still compatible within errors. 

The Q 2 dependence of the cross-section is shown in Figure |5[ There is a steep decrease 
in the cross-section with Q 2 . The expectation from the direct photon component only of the 
RAPGAP(LO)/DG model shows a rate of suppression of the cross-section which is independent 
of xif ts . However, the data show a rate of suppression which diminishes with increasing x^ ts . 
This behaviour is governed by the additional Q 2 suppression of the photon parton densities 
and the full RAPGAP(LO)/DG model including the resolved photon component gives a better 
description. 
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From the above comparisons we conclude that the observed dependence of the dijet cross- 
section is consistent with that predicted for a resolved virtual photon with parton density func- 
tions evolving with Q 2 according to QCD motivated models. In the next section we therefore 
proceed to extract an effective parton density for virtual photons from the data. 



6 The effective parton density for virtual photons 

6.1 Measurement of the Effective Parton Density 

In order to measure the parton densities of the virtual photons we adapt the Single Effective 
Subprocess Approximation [[3~T|], originally developed for use in pp collisions and recently used 
to investigate real photon structure []3p. In LO, the cross-section for the production of dijets by 
resolved virtual photons can be written as: 



d 5 



a 



dy dx y dx p d cos 6* dQ 2 

3 2 7T S on U Xrv Xr* 



£ Jlf/ W , ^ W ~7^t,«C , JWyt) |M.. (cosr) |2 . (15) 
k=T,L y ij 



Here fT, , and /j/ p are the densities of parton species % in transverse photons, longitudinal 
photons and the proton respectively. They are evaluated at the factorisation scale which we set 
equal to the renormalisation scale and choose to be P 2 . The are matrix elements for 2^2 
parton-parton hard scattering processes. The quantity s cp is the square of the centre of mass 
energy in the ep collision, 9* is the polar angle of the outgoing partons in the parton-parton 
centre of mass frame and x p is the momentum fraction of the parton out of the proton. The 
fluxes of transverse and longitudinal photons are given by [|3"2"1]: 



a 2(1 -y) 



a rl + (1 - y) 2 1 2m 2 e y 
2^ 



V Q 2 Q 4 



fj/e(y,Q ) = (16) 

The measured dijet cross-sections do not allow us to separate the various parton sub-processes 
or photon polarisation states. To obtain a factorisable form, we first approximate equation < JT3| ) 
by: 



d 5 



a 



dy dx 1 dx p d cos 6* dQ 2 

i ff /c (y,Q 2 ) f lh (x„p t 2 ,Q 2 ) f j/p (x p ,p t 2 



y J'W' 1 * ^ > \ Mi A C os6*)\ 2 , (17) 



327rs e p y . . x 7 Xp 

ij 

We have defined a set of photon polarisation-averaged parton densities: 

8 y 7 (x 7 , P 2 , Q 2 ) + e/ t y 7 



f lh (x„P 2 ,Q 2 ) = tf 7 (x 7 ,P 2 ,Q 2 ) + e/^ 7 (x 7 ,P 2 ,Q 2 ), (18) 
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where e ~ 1 is the ratio of longitudinal to transverse photon fluxes, averaged over the y-range 
of the data, fh is expected to be small over most of the kinematic range considered here [|7|, 
HI" 



The Single Effective Subprocess (SES) approximation exploits the fact that the dominant 
contributions to the cross-section, namely qg — > qg, gg — > gg and qq — > qq t-channel pro- 
cesses, come from parton-parton scattering matrix elements that have similar shapes and so 
differ mainly by their associated colour factors. Thus the sum over processes can be replaced 
by a single effective sub-process cross-section and effective parton densities for the photon and 
proton: 

d 5 a 



dy dx 7 dx p d cos 8* dQ 2 



i HMQ 2 ) U^PIQ*) UM1 lMms(cosrf . (19) 



327rs e p y Xp 
where the effective parton densities are 

; 7 (x 7) p 2 , q 2 ) = {M x v p l Q 2 ) + hh^ p l Q 2 )) + z M x v p h Q 2 ) C 20 ) 



and fp(x p ,P?) ee ^( /q/p (x p ,P 4 2 ) + /q /p (a: p ,P t 2 )) + ^/ g/p (x p) P 4 2 ) (21) 

and the sums are over the quark flavours. 

To extract the effective parton densities, the Bayesian unfolding method was applied again 
to correct the dijet cross-section to the diparton cross-section. This second unfolding corrects 
for hadronisation effects, the influence of the soft underlying event, and initial and final state 
QCD radiation. It uses correlations between the {x 7 e * s , E T } of the jets and {s 7 , P 2 } of the 
parton-parton hard scattering, obtained using the HERWIG(LO)-DG simulation. Here we use 
LO parton densities and the 1-loop formula for as for a consistent leading-order treatment. The 
resolution in x 7 varies from 15% at low x 7 to 10% at high x 7 in resolved events, and is 7% in 
direct processes. The P 2 resolution is 40% at low P t improving to 24% at high P t . The effective 
parton density was then determined by comparing the diparton cross-section measured in the 
data with that predicted by the simulations with a known set of photon parton densities: 

idata _ fMC x (d 3 a e p/dQ 2 dP 2 dx 7 ) DATA 

h h (dV ep /dQ 2 dP 2 dx 7 ) MC ; 



using equation ^0| to evaluate f^f c 



7 

Now we discuss the systematic errors on the effective parton density. Additional system- 
atic errors associated with the second unfolding were added in quadrature to those associated 
with determination of the triple-differential cross-section measurement described in the previ- 
ous section. The most significant new systematic effects arise from the model dependences. 
These were estimated by repeating the second unfolding with RAPGAP and HERWIG varying 
the input parton densities, the amount of soft underlying event and the hadronisation models as 
detailed below. 
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• Model Dependence. Additional model dependences arising from the use of different par- 
ton showering and hadronisation mechanisms and different input parton densities were 
estimated by unfolding the data using the RAPGAP(L0)/SAS-1D, RAPGAP(LO)/SAS- 
2D and RAPGAP(LO)/DG models and comparing the results with those obtained after 
unfolding with HERWIG(LO)/DG. We assign a further 25% systematic error on the basis 
of this test, of which 20% arises from the hadronisation uncertainties. 

• Unfolding Instability. Unfolding instabilities were estimated by varying the number of 
iterations and lead to a 10% uncertainty. 

• Soft Underlying Event. The transformation from jet-based observables to parton variables 
is more strongly influenced by the presence of the soft underlying event than was the case 
for the transformation between true and measured jets. As an independent measure of 
the amount of soft underlying event we examine the transverse energy flow in the region 
outside the jets. In Figure || we show the transverse energy flow per unit area in the (77, 4>) 
plane outside of circles of radius 1 .3 about the two highest E t jets and in the range — 1 .0 < 
77 < 1.0 as a function of Q 2 . This central region of pseudo-rapidity is where transverse 
energy flow from remnant-remnant interactions is expected to be largest. The energy flow 
is corrected for detector effects. The inner error bars show the statistical errors and the 
outer error bars the quadratic sum of statistical and systematic errors. The latter include all 
the sources considered for the jet cross-section measurement. The observed decrease with 
increasing Q 2 is compared with predictions from HERWIG(LO)/DG with 0%, 5% and 
10% soft underlying event and with RAPGAP(LO)/SAS-2D. Note that with the different 
choice of minimum P t and space-like shower parameter QSPAC (see section 3), which 
is used in HERWIG(LO)-DG for this step in the analysis, a lower soft underlying event 
frequency is required to describe the data. 

We varied the soft-underlying event probability in HERWIG between 0% and 10% in 
order to estimate the systematic uncertainty on the effective parton density. The resulting 
systematic error has a mean value of 20% and the error is largest at low x 7 and low Q 2 . 

6.2 Discussion of the Extracted Effective Parton Density 

The extracted effective parton density is given in Table |] and shown in Figures [7], H and |9|. The 
measurement uncertainties are everywhere domininated by the systematic errors. Only points 
where (P 2 ) > (Q 2 ) are shown. It should be noted that in the highest Q 2 bins, Q 2 approaches 
P 2 and the assumptions of factorisation involved in the definition of a universal effective parton 
density begin to break down [0, [28p. Furthermore, the effects arising from higher twist con- 
tributions and longitudinally polarised photons are largest here. Nevertheless, we extract the 
effective parton distributions in the full region of the available data. The universality of the 
PDF's extracted in the region where Q 2 is of the same order as P 2 will need to be demonstrated 
in other virtual photon induced reactions. The Figures show the evolution of the effective parton 
density in the photon both with the scale at which it is probed and with its virtuality. The data 
are compared with three sets of parton densities, SAS-1D, SAS-2D and the effective parton 
density calculated from the Drees-Godbole model using GRV-LO densities for the real photon 
and setting the free parameter, to, in the suppression factor to 0.1 GeV. The predictions were 
evaluated at the mean x 7 and logarithmic mean Q 2 and P 2 of the ranges. 
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Figure [7] shows that the effective parton density tends to rise with x 7 in the region studied. 
This shape is characteristic of photon structure. The data are described by all three models 
within errors except possibly in the two higher Q 2 ranges where the models tend to underesti- 
mate the data in the intermediate and high s 7 bins. 

In Figure ^ the parton density is shown as a function of the square of the probing scale P 2 
in ranges of x 7 and Q 2 . The effective parton density is roughly independent of P 2 in each of the 
ranges. This scale dependence contrasts with the falling behaviour expected for hadronic parton 
densities and, within the rather large systematic errors, is consistent with the normalisation and 
logarithmic scale dependence predicted by perturbative QCD for the anomalous component of 
the photon. 

The decrease of the parton density with virtuality is most clearly seen in Fig |9|, where the 
Q 2 dependence is shown in ranges of P 2 and x 1 . The three parameterisations for the parton 
density all give a good description of the data both in the lowest x 7 range and in the lowest two 
Q 2 bins. They all predict a more rapid suppression as Q 2 increases than is seen in the data. It 
is in this region where Q 2 — ► P 2 that non-leading terms, not accounted for in the models, are 
expected to become important and may affect the extraction of the effective parton distribution 
from these data. 

In Figure [10|, the Q 2 evolution of the effective parton density is shown for P 2 = 85 GeV 2 , 
in the upper Figure for 0.35 < x 7 < 0.5 and in the lower for 0.5 < x 7 < 0.7. Superimposed 
on the same plots are photoproduction data at P 2 = 112 GeV 2 and x 7 values of 0.3 and 0.55, 
taken from reference [TJ], which we have extrapolated to the P 2 and x 7 values in the Figures. 
The extrapolation was based on the GRV-LO parton densities which give a good description of 
the photoproduction data. The evolution is compared with that predicted by the three models 
described above and also with a simple p-pole suppression factor characteristic of a pure VMD 
model: 

/,| 7 (x 7 ,P 2 ,Q 2 ) ~ / i|7 (x 7 ,P 2 ,0)(-^^) 2 . (23) 

The latter clearly underestimates the data. The logarithmic suppression predicted by the virtual 
photon models on the other hand gives a good description below Q 2 ~ 20 GeV 2 . At higher Q 2 
they predict a more rapid decrease than is seen in the data. 



7 Conclusions 

We have measured the triple-differential cross-section, d 3 a ep /dQ 2 dE t dxlf^, in dijet events 
for 1.6 < Q 2 < 80 GeV 2 and 0.1 < y < 0.7. The measured cross-sections show the x^ ts , 

E t and Q 2 behaviour expected for processes in which a virtual photon, carrying a partonic 
structure evolving according to perturbative QCD, interacts with the proton via hard parton- 
parton scattering. The measurements are consistent with the perturbative QCD prediction that, 

2 

as Q 2 — > E t , the photon structure reduces to a simple direct coupling to qq pairs and the dijet 
cross-section is well-described without invoking photon structure. LO Monte Carlo models 
based on these QCD predictions give a good description of the data. 
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An effective, leading order, parton density of the virtual photon has been extracted in the 
Single Effective Subprocess approximation and its dependences on x 7 , probing scale P 2 and 
target virtuality Q 2 have been measured. The photon parton density is approximately indepen- 
dent of and, within errors, it is consistent with the normalisation and logarithmic scaling 
violations characteristic of photon structure. It is seen to be suppressed with increasing Q 2 as 
predicted by perturbative QCD. 
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Table 2: The triple-differential dijet cross-section d 3 a ep /dQ 2 dE t dx^ ets for 0.1 < y < 0.7 in 
ranges of Q 2 , E t and x^ ts . The cross-section in pb is given together with the statistical, positive 
systematic and negative systematic errors. The table shows measurements for Q 2 < 8.0 GeV 2 . 
The higher Q 2 region is given in table ^. 
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Table 3: The triple-differential dijet cross-section d 3 <7 ep /dg 2 d-E t dx 7 ets in ranges of Q 2 , E t 
and xi^ ts . The cross-section in pb is given together with the statistical, positive systematic and 
negative systematic errors. The table shows measurements for Q 2 > 8.0 GeV 2 . The lower Q 2 
region is given in table |[ 
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Table 4: The leading order effective parton density of the photon x 1 f 1 = ^2 (/ q / 7 + /q/ 7 ) + 
| / g / 7 , divided by the fine structure constant a, for different values of Q 2 , P 2 and x 7 . Statistical, 
positive systematic and negative systematic errors are given. 



22 




Figure 1 : Examples of leading order Feynman diagrams of dijet production in 7*p collisions. In 
direct processes (left) all of the photon's 4-momentum enters the hard collision. If the transverse 
momentum of the outgoing dijets (labelled jet 1 and jet 2 in the figures) is large compared with 
Q 2 , the virtual photon may be resolved (right) with only a fraction x y < 1 of the photon's 4- 
momentum involved. The diagrams are applicable to dijet production with either real or virtual 
photons. 



23 



H1 Data 



-3.5<?}<-2.0 



HERWIG(HO)/DG, 10% SUE 
— RAPGAP(HO)/SAS-2D 

-2.0<tj<-1 .25 -1.25<7j<0.0 



- 

h» * 

i i i i i 


• 

• 

• 

• 

i i i i i 


Till 


• 

* 

1 1 1 1 1 


: J 

iii 


i 

• '» 

iii 


- 

I* • 

i |» 

: ^ C» • 

iL *•[ 

i i i i i 


* 

* 

• 

• 

i i i i i 


I n 
.. p 

_ L=j *WW 

1 1 1 1 1 


1 

• 

• 

i i i i i 


[ J 

I* « 

i i i i i 


l 

— 1 

i i i i i 


•- 

« 

* r 
J* 

i i -P -1 1 i 


• 

• 

• 

• 

r 

iii 


f* 

In * 
* 

Fi * 
— I. J* r 

- 

iii 


i 

• 

iii 


-\ r 

r I* 
1* 

h» » 

iii 


> 

•i 

• 



-2.5 



2.5 -2.5 



2.5 -2.5 



2.5 



Figure 2: The observed transverse energy flow in with respect to the axes (located at = 0) 
of each of the two highest E t jets in the selected event sample for a slice | rj — rf et |< 1. The 
result is shown for various ranges of 77 and E t of the jets. The second jet is always chosen to 
have 4> < 0, leaving the pedestal level clearly visible to the right of the jet core, in the region 
> 0. The data are compared with predictions from the HERWIG(HO)/DG simulation with 
10% soft underlying event (dashed histogram) and with the RAPGAP(HO)/SAS-2D model. 
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---- HERWIG(LO)/DG cj = 0.2 
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Figure 3: The differential dijet cross-section d 3 a ep /dQ 2 dE t dxi^ ts shown as a function of x^ ts 

for different regions of E t and Q 2 . The scale factors applied to the cross-sections are indicated. 
The error bar shows the quadratic sum of systematic and statistical errors. The absence of a 
data point indicates that no measurement was made because of insufficient statistics for the two 
dimensional unfolding. Also shown is the HERWIG(LO)/DG model with 10% soft-underlying 
event and two choices of the Q 2 suppression factor cu. The full histogram is for uj = 0.1 GeV 
and the dashed for cu = 0.2 GeV. The direct component of this model is shown as the shaded 
histogram. 
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E r /GeW 2 



Figure 4: The differential dijet cross-section d 3 a ep /dQ 2 dE t dxi^ ts shown as a function of E t 
for different regions of x^ ts and Q 2 . The error bars show the quadratic sum of systematic and 
statistical errors. Also shown is the RAPGAP(LO) model with three choices of photon parton 
density. The DG model with GRV-LO real photon densities and uo = 0.2 GeV is shown as the 
full curve. The predictions with SAS-1D and SAS-2D are shown as the dashed and dot-dashed 
curves respectively. 
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RAPGAP(LO)/DIR • H1 Data 
RAPGAP(LO)/DG 

0.15 < xf < 0.3 0.3 < x T ie,s < 0.45 0.45 < x y ie,s < 0.6 0.6 < x r icts < 0.75 0.75 < x r iets < 1.0 




1 10 



Q 2 /GeV 2 

Figure 5: The differential dijet cross-section d 3 a ep /dQ 2 dE t dx 3 7 ets shown as a function of Q 2 

for different regions of x 3 y ets and E t . The error bars show the quadratic sum of systematic and 
statistical errors. The prediction from the RAPGAP(LO)/DG model is shown as the full curve. 
The dashed curve shows the prediction for the direct photon processes in this model. 
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Figure 6: The average transverse energy per unit area in rj— <fi space in the range —1.0 < rj < 1.0 
and < < 2n in the 7*p cms as a function of Q 2 and outside the two highest E t jets. The 
data are corrected for detector effects. The inner error bars show the statistical errors and the 
outer error bars are the quadratic sum of the statistical and systematic errors. Also shown 
is the prediction from the HERWIG(LO)/DG model with three choices for the percentage of 
resolved events with a soft underlying event: 0% (dash-dotted), 5% (full) and 10% (dashed). 
The prediction from the RAPGAP(LO)/SAS-2D model is shown as the dotted curve. 
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Figure 7: The leading order effective parton density of the photon x 1 f 1 = ^ (/ q / 7 + /q/ 7 ) + 
| / g / 7 , divided by the fine structure constant a, as a function of x 7 for different values of Q 2 
and P t 2 . The data are displayed as points, with the inner error bar depicting the statistical error, 
and the total error bar the quadratic sum of statistical and systematic errors. In most bins the 
inner errors are contained within the data point marker. Also shown are the predictions from 
the DG model using GRV-LO real photon parton densities and u — 0.1 GeV (solid line) and 
the SAS-1D (dashed line) and SAS-2D (dot-dashed line) parameterisations. 
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Figure 8: The leading order effective parton density of the photon x 7 / 7 = ^ (/ q / 7 + /q/ 7 ) + 
| / g / 7 , divided by the fine structure constant a, as a function of the squared parton transverse 
momentum, P t 2 , for different values of Q 2 and x 7 . The data are displayed as points, with 
the inner error bar depicting the statistical error, and the total error bar the quadratic sum of 
statistical and systematic errors. In most bins the inner errors are contained within the data 
point marker. Also shown are the predictions from the DG model using GRV-LO real photon 
parton densities and uo = 0.1 GeV (solid line) and the SAS-1D (dashed line) and SAS-2D 
(dot-dashed line) parameterisations. 
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Figure 9: The leading order effective parton density of the photon x 1 f 1 = ^ (/ q / 7 + /q/ 7 ) + 
| / g / 7 , divided by the fine structure constant a, as a function of Q 2 for different values of Pf 
and x 7 . The data are displayed as points, with the inner error bar depicting the statistical error, 
and the total error bar the quadratic sum of statistical and systematic errors. In most bins the 
inner errors are contained within the data point marker. Also shown are the predictions from 
the DG model using GRV-LO real photon parton densities and u — 0.1 GeV (solid line) and 
the SAS-1D (dashed line) and SAS-2D (dot-dashed line) parameterisations. 
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Figure 10: The leading order effective parton density of the photon x 1 f 1 = (/ q / 7 + /q/ 7 ) + 
| / g / 7 , divided by the fine structure constant a, as a function of Q 2 for P t 2 = 85 GeV 2 and two 
values of x 1 . The data are displayed as points, with the inner error bar depicting the statistical 
error, and the total error bar the quadratic sum of statistical and systematic errors. The points 
at Q 2 = are taken from reference [0j and extrapolated to the right P 2 and x 7 values by 
scaling with factors derived from GRV-LO parton densities for real photons. Also shown are 
the prediction from the DG model using GRV-LO photon parton densities and lo = 0.1 GeV 
(solid line) and the SAS-1D (dashed line) and SAS-2D (dotted line) parameterisations. The 
dot-dashed curve shows the photoproduction data scaled by a p-pole factor (see text). 
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